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ABSTRACT 

We quantitatively study the probability distribution function (PDF) of cos- 
mological nonlinear density fluctuations from N-body simulations with Gaussian 
initial condition. In particular, we examine the validity and limitations of one- 
point and two-point log-normal PDF models against those directly estimated 
from the simulations. We find that the one-point log-normal PDF describes very 
accurately the cosmological density distribution even in the nonlinear regime (the 
rms variance a n \ ^ 4 and the over-density 5 ^ 100). Furthermore the two-point 
log-normal PDFs are also in good agreement with the simulation data from linear 
to fairly nonlinear regime, while slightly deviate from them for 5 ^ — 0.5. Thus 
the log-normal PDF can be used as a useful empirical model for the cosmological 
density fluctuations. While this conclusion is fairly insensitive to the shape of the 
underlying power spectrum of density fluctuations P{k), models with substantial 
power on large scales, i.e., n = d In P(k)/d In k £ — 1, are better described by the 
log-normal PDF. On the other hand, we note that the one-to-one mapping of the 
initial and the evolved density fields consistent with the log-normal model does 
not approximate the broad distribution of their mutual correlation even on av- 
erage. Thus the origin of the phenomenological log-normal PDF approximation 
still remains to be understood. 



Subject headings: cosmology: theory - galaxiesxlustering - galaxies: dark matter 
- large-scale structure of universe - methods: numerical 
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1. INTRODUCTION 

Probability distribution function (PDF) of the cosmological density fluctuations is the 
most fundamental statistic characterizing the large-scale structure of the universe. In the 
standard picture of gravitational instability, the PDF of the primordial density fluctuations 
which are responsible for the current structures in the universe is assumed to obey the 
random-Gaussian. Therefore it is fully specified by the two-point correlation function £(r), 
or equivalently, the power spectrum P{k). As long as the density fluctuations are in the 
linear regime, their PDF remains Gaussian. Once they reach the nonlinear stage, however, 
their PDF significantly deviates from the initial Gaussian shape due to the strong non- 
linear mode-coupling and the non-locality of the gravitational dynamics. The functional 
form for the resulting PDFs in nonlinear regimes are not known exactly, and a variety of 
phenomenological models have been proposed (Saslaw 1985; Suto, Itoh, & Inagaki 1990; 
Lahav et al. 1993; Gaztahaga & Yokoyama 1993; Suto 1993; Ueda & Yokoyama 1996). 
Once such one-point PDF is specified, one can characterize the clustering of the universe 
with the higher-order statistics like skewness and kurtosis. Moreover the two-point PDF is 
useful in estimating the errors in the one-point statistics due to the finite sampling since the 
measurement at different positions is not independent and their correlations are supposed to 
be dominated by the two-point correlation function (Colombi, Bouchet, & Schaeffer 1995; 
Szapudi & Colombi 1996). Also the two-point PDF plays an important role in analytical 
modeling of dark halo biasing on two-point statistics. 

From an empirical point of view, Hubble (1934) first noted that the galaxy distribution 
in angular cells on the celestial sphere may be approximated by a log-normal distribution, 
rather than a Gaussian. More recent analysis of the three-dimensional distribution of galaxies 
indeed confirmed this (e.g., Hamilton 1985; Bouchet et al. 1993; Kofman et al. 1994). 
Interestingly, several N-body simulations in cold dark matter (CDM) models also indicated 
that the PDF of density fluctuations is fairly well approximated by the log-normal (e.g., Coles 
& Jones 1991; Kofman et al. 1994; Taylor & Watts 2000), at least in a weakly nonlinear 
regime. 

Those observational and numerical indications have not yet been understood theoreti- 
cally; Bernardeau (1992, 1994) showed that the PDF computed from the perturbation theory 
in a weakly nonlinear regime approaches the log-normal form only when the primordial power 
spectrum is proportional to k n with n = — 1. On the basis of this result, Bernardeau & Kof- 
man (1995) argued that the successful fit of the log-normal PDF in the CDM models should 
be interpreted as accidental, and simply resulted from the fact that those model have the 
power spectrum well approximated by k UcB with n e fr — — 1 on scales of cosmological interest. 
They claimed that the log-normal PDF may fail either in a highly nonlinear regime or in 
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models with power spectrum with n c s ^ — 1. In fact, Ueda & Yokoyama (1996) conclude 
that the log-normal PDF does not fit well the PDF in a highly nonlinear regime, from the 
analysis of CDM simulations by Suginohara & Suto (1991) employing N = 64 3 particles in 
a lOOMpc box. 

The aim of this paper is to study the extent to which the log-normal model describes the 
PDF in weakly and highly nonlinear regimes using the high-resolution N-body simulations 
with iV = 256 3 . In particular, we extend our analysis to the two-point PDF, in addition 
to the one-point PDF discussed previously. Bernardeau (1996) analytically computed the 
two-point PDF using the perturbation technique and compared somewhat indirectly with N- 
body simulations in a weakly nonlinear regime. In contrast, we focus on the highly nonlinear 
regime, and examine the validity of the empirical log-normal model. 

This paper is organized as follows. Section 2 describes the log-normal PDF derived 
through the one-to-one mapping between the linear and nonlinear density fields. The detailed 
comparison between the log-normal predictions and N-body results is presented in §3. Finally 
§4 is devoted to conclusions and discussion. 



2. PROBABILITY DISTRIBUTION FUNCTIONS FROM THE 
LOG-NORMAL TRANSFORMATION 

In this section we briefly outline the derivation of the log-normal PDF assuming the 
one-to-one corresponding between the linear and evolved density fluctuations. Throughout 
the paper, we consider the mass density field, p(x; R) at the position x smoothed over the 
scale R. This is related to the unsmoothed density field p(x) as 

p{x-R) = J d 3 yW(\x-y\;R) p(y) 

= J ^W{kR)p{k)e'^ x . (1) 



In the above expression, W denotes the window function, and W and p represent the Fourier 
transforms of the corresponding quantities. In what follows, we adopt the two conventional 
windows: 

[3(sini-xcosxj/r (Top-hat). 

Then the density contrast at the position x is defined as 8(x; R) = (p(x; R) — p)/p, with p 
denote the spatial average of the smoothed mass density field. For simplicity we use 5 to 
denote 5(x; R) unless otherwise stated. 
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2.1. One-point Log-normal PDF 

The one-point log-normal PDF of a field 5 is defined as 

{ln(f + 5) + a 2 /2} 2 



yj2lt<j( 



2o\ 



1 + 6 



(3) 



The above function is characterized by a single parameter a x which is related to the variance 
of 5. Since we use 5 to represent the density fluctuation field smoothed over R, its variance 
is computed from its power spectrum P Q \ explicitly as 

1 f 00 

a 2 nl (R) = —j^ P IlX {k)W 2 {kR)k 2 dk. (4) 

Here and in what follows, we use subscripts "lin" and "nl" to distinguish the variables corre- 
sponding to the primordial (linear) and the evolved (nonlinear) density fields, respectively. 
Then <j\ depends on the smoothing scale R alone and is given by 

a 1 2 ( J R) = ln[l+a r 2 ll (i?)]. (5) 

Given a set of cosmo logical parameters, one can compute <7 n i(i?) and thus <J\{R) very accu- 
rately using a fitting formula for P n \(k) (e.g., Peacock & Dodds 1996, hereafter PD). In this 
sense, the above log-normal PDF is completely specified and allows the definite comparison 
against the numerical simulations (§3). 

It is known that the above log-normal function may be obtained from the one-to-one 
mapping between the linear random-Gaussian and the nonlinear density fields (e.g., Coles 
& Jones 1991). Define a linear density field g smoothed over R obeying the Gaussian PDF: 

411(9) ^vir p (-4)' (6) 

where the variance is computed from its linear power spectrum: 

al(R) = — ^ PUk)W 2 (kR)k 2 dk. (7) 

If one introduces a new field S from g as 

1 f 9 



1 + 6= -== exp f- + a 2 J , (8) 



the PDF for 5 is simply given by {dg / dS)PQ\g) which reduces to equation (3). 
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At this point, the transformation (8) is nothing but a mathematical procedure to relate 
the Gaussian and the log-normal functions. Thus there is no physical reason to believe 
that the new field 5 should be regarded as a nonlinear density field evolved from g even in 
an approximate sense. In fact it is physically unacceptable since the relation, if taken at 
face value, implies that the nonlinear density field is completely determined by its linear 
counterpart locally. We know, on the other hand, that the nonlinear gravitational evolution 
of cosmological density fluctuations proceeds in a quite nonlocal manner, and is sensitive to 
the surrounding mass distribution. 

Nevertheless the fact that the log-normal PDF provides a good fit to the simulation 
data empirically as discussed in §1 implies that the transformation (8) somehow captures an 
important aspect of the nonlinear evolution in the real universe. In §3, we present detailed 
discussion on this problem. Before that, we derive the two-point log-normal PDF by applying 
this transformation in the next subsection. 



2.2. Two-point Log-normal PDF 

Consider two density fields, 5\ = S(xi,R) and S 2 = 5(x2',R), located at X\ and x 2 
smoothed over R. We denote the two-point PDF by P^ 2 \5i, 5 2 ; r) for the two fields with a 
specified separation r, i.e., satisfying the condition \x 1 — x 2 \ — r. 

In the case of the Gaussian fields gi and g 2 , this two-point PDF is given by the bi-variate 
Gaussian (e.g., Bardeen et al. 1986): 



(9) 



where 



and 



V (9192) (gl) J V 6m(r) < ; 

6in(r ;j R) = ^ jT ' P lh ^k)W\kR)*-^^k 2 dk. (11) 

From an analogy of equation (8), let us assume that the transformation from (gi,g 2 ) to 
(Si, S 2 ) is given by the form: 

l + 5 i = ae l3g \ (i = l,2) (12) 
The coefficients a and (3 are determined by the following conditions: 

(*i> = <*2> = 0, (13) 
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<*?> = (S 2 2 ) 



J2 



(14) 

(Si5 2 )=U(r). (15) 

In the above expressions, the angular bracket denotes the average over the two-point PDF 
which in the present model, reduces to 

(F(5 U 5 2 )) = F(6 1 ,6 2 )pW(6 1 ,6 2 ;r)d6 1 d5 2 

p poo 

= // ^{^i{9i)^2(g2))PG(gi,92;r)dg 1 dg 2 . (16) 

J J — oo 

After a straightforward calculation, one obtains 

/m(l 



a 



slin 



Then this procedure yields the two-point log-normal PDF: 



P^(5 1 ,5 2 ;r) = 



1 



2nVS 2 - X 2 



exp 



S{L\ + L\] — 2XLiL 2 



2{S 2 - X 2 } 



(l + 5 1 )(l + 5 2 y 



where 



X EE ln(l + £nl), 

5 ee ln(l + 4), 



Li = In |(l + ^)y / l+^|, (i = l,2). 



Again the nonlinear two-point correlation function C, Q \(r; R) can be computed as 



2ir 2 



Pni(k)W 2 {kR) 



sin (Ax) 
kr 



k 2 dk, 



(17) 

(18) 

(19) 
(20) 

(21) 
(22) 



and thus equation (18) can be fully specified using the PD nonlinear power spectrum. 



3. THE LOG-NORMAL PDFS AGAINST N-BODY SIMULATIONS 

The previous section discussed a prescription to derive one-point and two-point log- 
normal PDFs assuming the one-to-one mapping between the linear and the nonlinear density 
fields. As remarked, however, the assumption does not seem to be justified in reality. So in 
this section we compare the log-normal PDFs extensively with the results of cosmological N- 
body simulations, and discuss their validity and limitations. The analysis for the one-point 
PDF below significantly increases the range of S compared with several previous work. As 
far as we know, the direct estimation of the two-point PDF in the nonlinear regime from 
simulations has not been performed before and this is the first attempt. 
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3.1. N-body Simulations 

For the present analysis, we use a series of cosmological N-body simulations in three 
CDM models (SCDM, LCDM and OCDM for Standard, Lambda and Open CDM models, 
respectively; Jing & Suto 1998) and four scale-free models with the initial power spectrum 
P(k) oc k n (n=l, 0, -1, and -2; Jing 1998). All the models employ N = 256 3 dark 
matter particles in a periodic comoving cube Ll ox , an d are evolved using the P 3 M code. The 
gravitational softening length is L box /2560 (3L box /5120) for the CDM (scale-free) models, 
and kept fixed in the comoving length. The amplitude of the fluctuations in CDM model , 
cr 8 , is normalized according to the cluster abundance (e.g., Kitayama & Suto, 1997). The 
scale- free models assume the Einstein-de Sitter universe (the density parameter flo — l, the 
cosmological constant Ao = 0). The other parameters of the CDM models are summarized 
in Table 1. 

The mass density fields are computed on 512 3 grids in the simulation box. First we 
assign particles to each grid point using the cloud-in-cell interpolation. Then we apply the 
smoothing kernel in the Fourier space, and then obtain the smoothed density fields after the 
inverse Fourier transform. Note that the density fields on those grids are not completely 
independent for R > L box /512, and we do heavy over-sampling in this sense. Nevertheless 
the error-bars quoted in our results below are estimated from the variance among the three 
different realizations for each model (except the n — — 1 model which has two realizations 
only), and thus are free from the over-sampling. 

3.2. The one-point PDF 

Consider first the one-point PDFs in CDM models (Fig. 1). The PDFs are constructed 
by binning the data with AS = 0.1, but we do not plot all the data points just for an 
illustrative purpose. We compute the density fields smoothed over Gaussian (Left panels) and 
Top-hat (Right panels) windows with different smoothing lengths;/? = 2/i -1 Mpc, 6/i _1 Mpc 
and 18/i _1 Mpc plotted in cyan, red and green symbols with error-bars, respectively. The 
corresponding values of o~ n \ are summarized in Table 2, and also shown on each panel. Solid 
lines show the log-normal PDFs adopting those a n \ directly evaluated from simulations. The 
agreement between the log-normal model and the simulation results is quite impressive. A 
small deviation is noticeable only for 5 £ — 0.5. 

We also show the log-normal PDFs in dashed lines, adopting cr n i calculated from the 
nonlinear fitting formula of PD (values in parenthesis in Table 2). Therefore the predictions 
do not use the specific information of the current simulations, and are completely indepen- 
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dent in this sense. While these predictions are in good agreement with simulation data for 
R £j 6/i -1 Mpc, the results for R = 18/i _1 Mpc are rather different. Actually this discrep- 
ancy should be ascribed to the simulations themselves, not to the model predictions; Table 
2 indicates that the a n \ in the current simulations become systematically smaller than the 
PD predictions for larger R. This is because the simulations assume (incorrectly) no fluc- 
tuations beyond the scale of the simulation boxsize Lbox- This constraints systematically 
reduce the fluctuations as the smoothing scale R approaches L hox . We made sure that this 
is indeed the case by repeating the same analysis using the CDM simulations evolved in 
^box = 300/i -1 Mpc (Jing & Suto 1998); the variance of fluctuations at R < 18/i _1 Mpc from 
the simulations agrees with the PD prediction within 2% accuracy. Thus we conclude that 
the log-normal PDF with the PD formula reproduces accurately the simulation results in 
the CDM models. 

Next turn to the scale-free models. Figure 2 shows the similar plots corresponding to 
Figure 1 but for n — 1 to —2 models (from top to bottom panels). In this figure, we compare 
the simulation data (symbols with error-bars) with the log-normal PDF predictions (solid 
lines) adopting a n \ from simulations (Table 3). Generally their agreement is good also in 
these models. A closer look at Figure 2, however, reveals that the simulation results start 
to deviate from the log-normal predictions at both high and low density regions, and that 
the deviation seems to systematically increase as n becomes larger. While this tendency 
is qualitatively consistent with the earlier claim by Bernardeau (1994) and Bernardeau & 
Kofman (1995) on the basis of the perturbation theory, our fully nonlinear simulations show 
that the deviation from the log-normal PDF is not so large even in these scale-free models. 

To examine the validity of the log-normal PDF more quantitatively, we compare the 
normalized skewness S = (S 3 )/(S 2 ) 2 and the normalized kurtosis K = ((5 4 ) — 3(5 2 ) 2 )/ (S 2 ) 3 . 
The log-normal PDF predicts that 



which are plotted in dotted lines in Figure3 for six models; two LCDM models with L box = 
100/i _1 Mpc and 300/i _1 Mpc, and four scale-free models with n — 1, 0, -1, and —2. 

In practice, however, the density field 5 in numerical simulations does not extend the 
entire range between —1 and oo, but rather is limited as <5 m i n < 5 < 5 max due to the finite 
size of the simulation box. Thus the n-th order moments of S in simulations may be better 
related to 



S(R) 
K{R) 



S + a 2 nl (R), 

16 + 15a^(R) + 6a^(R) + a^(R), 



(23) 
(24) 
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The specific values for <5 min and S max may be roughly estimated from the condition that the 
expectation number of independent sampling spheres in the simulation box for 5 < 5 min or 
8 > <5 max becomes unity: 



Dashed lines in Figure 3 show the log-normal PDF predictions based on equations (25) 
to (27). The filled triangles and squares represent the measurement of S and K from the 
simulations, and finally the solid lines indicate the log-normal PDF predictions using equation 
(25) with the actual values for <5 min and <5 max in the simulations. Except for the n = 1 scale- 
free model, the predictions in solid lines reproduce the simulation data very well, which 
indicates that the log-normal PDF is in fact a good approximation. The relatively large 
discrepancy between the log-normal prediction and the simulation in the n = 1 model is 
real since one can clearly recognize the systematic tendency with respect to n; models with 
smaller n, i.e., with substantial power on large scales , are better described by the log-normal 
PDF. This is consistent with the discussion by Bernardeau (1994). 

Incidentally both the current simulations and the log-normal PDF approximation con- 
firmed the relatively strong scale-dependence of S and K for a n \ > 1 as pointed out earlier 
by Lahav et al. (1993) and Suto (1993). In fact the degree of their scale-dependence is also 
sensitive to the underlying power spectrum. Thus the hierarchical ansatz for the higher-order 
clustering is not valid in general. 

In summary, we find that the one-point log-normal PDF remains a fairly accurate model 
for the cosmological density distribution even up to a n \ ~ 4 and 5 ~ 100, fairly independently 
of the shape of the underlying power spectrum of density fluctuations. The range of validity 
turns out to be significantly broader than those from the previous studies based on lower- 
resolution simulations; 0.1 0.6 and 5 < 4 (Kofman et al. 1994), and 0.3 1.5 and 
5^9 (Bernardeau & Kofman 1995), for instance. 



While we would like to perform the similar comparison for the two-point PDFs, it is 
a function of four variables, 5 1 , 5 2 , R, and r, and thus the comparison becomes rather 




(26) 



(27) 



3.3. The two-point PDF 
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complicated. Therefore we decided to use the conditional two-point PDF for the purpose: 

(28) 



1 11 2 ' } ~ P0){8 2 ) ' 



Since we already made sure that the one-point PDF is very accurately approximated by the 
log-normal, our task is to see if simulation results fit the conditional two-point log-normal 
PDF: 



S 



2ir(S 2 - X 2 ) 



cxp 



{SL 1 -XL 2 f 
' 2S(S 2 - X 2 ) 



1 



(29) 



The evaluation of the conditional two-point PDFs from simulations is carried out as 
follows. From the smoothed density fields computed on the 512 3 grid points, we first select 
those grid points with 5 2 — A5 2 /2 < 5 < 5 2 + A5 2 /2. The bin size AS 2 is adjusted for each 
value of 5 2 so that approximately 10 5 grid points satisfy the condition. Next we pick up 
all grid points separated at r — Ar/2 ~ r + Ar/2 from the above grids. The separation 
interval Ar is chosen so that 47rr 2 Ar(5i2/L b0 x) 3 ~ 10 3 - Finally we compute the conditional 
two-point PDF with a constant bin size of 0.1 in 8i. 

Figure 4 plots the resulting PDFs for the LCDM model with the Gaussian smoothing 
window; the upper four panels show the PDFs for the separation r = 4 and 6/i _1 Mpc and the 
smoothing length R = 2/i -1 Mpc, while the lower four panels for r = 12 and 18/i _1 Mpc and 
R = 6/i -1 Mpc. Solid lines indicate the conditional log-normal PDFs adopting the values of 
<j n i and £ n i from simulations, while dashed lines show those using the PD predictions (Tables 
2 to 5). Clearly the log-normal PDF is a reasonably good approximation. The deviation at 
^2 ~ — 0.5, on the other hand, seems real and may be an enhanced feature that we noted in 
the one-point PDF (Fig. 1). 

Figure 5 indicates that the good agreement is achieved not only in the Gaussian smoothed 
LCDM model, but also can be found in the other models and/or the top-hat smoothing. Fig- 
ure 6 plots the conditional two-point PDFs at 5 2 —0.7 and 5 2 ^ 10 in the scale- free models, 
where the deviation from the log- normal becomes manifest. Considering the error-bars esti- 
mated from the different realizations for each model, the deviation seems statistically real. 

As in the case of the one-point PDF, we illustrate the validity of the two-point log- 
normal PDF using the moments. Specifically we evaluate ((^i^) 2 ) and ((^i^) 3 ) according 
to 

((W)(r)= II (S 1 S 2 ) n P^(S 1 ,S 2 ;r)dS 1 dS 2: (30) 

JJ C(S!,S 2 ) 
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where we select the range of the integration as 

C(5 1 ,5 2 ) = {(Si,S 2 )\S min < 5i < 5 max ,5 min < 5 2 < 5 max }, (31) 

from the values of £ min and <5 max directly measured from each simulation model. The results 
are summarized in Table 6, which indicates again that the two-point log-normal PDF pre- 
dictions reproduce the simulation data except for n > 0. Thus we also conclude that the 
two-point log-normal model describes fairly well the PDF of cosmological fluctuations for 
most regions of 5i and 5 2 of interest; the small but finite deviations exist only in 5i ?i — 0.5 
and/or 5 2 ~ — 0.5, and in 5i ^ 10 and 5 2 ^ 10. 



3.4. Does the log-normal transformation approximate the gravitational 

evolution of the density fluctuations ? 

The agreement between the log-normal predictions and the simulation results might 
be interpreted as indirect evidence that the log- normal transformation (eq. [8]) is a good 
approximation for the nonlinear gravitational growth of the cosmological density fluctuations, 
at least on average. 

In order to see if this is really the case, we consider the relation of the smoothed density 
fields at the same comoving position but at different redshifts z\ and z 2 . For this purpose, we 
use one realization from the LCDM model evolved in L^ ox = 300/i _1 Mpc. If the log-normal 
transformation (8) is exact, the density fluctuations, 5(z\) and S(z 2 ), should satisfy 



1 + 6 ^ = h\ ( \ exp 
V 1 + W2) 



VMi + wi)} l v 



(32) 



Figure 7 plots the color contour of the joint probability P(8(z),8(z = 9)) of densities 
at z — 2.2, 1.0, and against that at z = 9 on the same grid points in the LCDM model. 
We adopt the Gaussian smoothing with R = 6/i -1 Mpc (Left) and 2/i~ 1 Mpc (Right). The 
solid lines in white and magenta represent the the log-normal transformation (32) and the 
conditional mean from simulations for a fixed 5(z = 9). The log- normal transformation 
traces the mean relation of simulations to some extent only when the nonlinearity is weak 
(see, higher z and larger R cases). On the other hand, in the nonlinear region, the trans- 
formation (32) starts to deviate from the mean relation of the simulations significantly, and 
the distribution around the mean relation becomes broad. The similar tendency was found 
in a somewhat different analysis by Coles, Melott, & Shandarin (1993). In a sense this is a 
physically natural and expected result, but then it makes even more difficult to account for 
the good agreement between the log-normal and simulation PDFs in those scales. 
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To better understand the distribution of the linear and its evolved density fields, we 
compute the conditional probability P(S(z) \ 5(z = 9)), i.e., the slice of Figure 7 at a 
given 5(z = 9). The results are plotted in Figure 8 which exhibits a some regularity in the 
distribution. The peak positions seem to show some scaling with respect to the value of 
5(z = 9), and also the tail of the distribution asymptotically approaches a single power-law. 
While we do not yet fully understand the behavior, this regularity in the distribution function 
may be useful in explaining our findings that the one-point and two-point log-normal PDFs 
work well empirically. 



4. CONCLUSIONS AND DISCUSSION 

In the present paper, we have estimated the probability distribution functions of cosmo- 
logical density fluctuations from the high-resolution N-body simulations with the Gaussian 
initial condition. In particular, we have critically examined the validity of the log-normal 
models for the one- and two-point PDFs both in weakly and in strongly nonlinear regimes. 

We have shown that the one-point log-normal PDF is a fairly accurate model not only 
in a weakly nonlinear regimes as claimed previously, but also in more strongly nonlinear 
regimes even up to a n \ ~ 4 and 5 ~ 100. Furthermore, we extended the analysis to the 
two-point PDF, and found that the log-normal PDF serves also as an empirically accurate 
model for the range of densities of interest. This is the case fairly independently of the shape 
of the underlying power spectrum of density fluctuations, although models with large power 
on small scales (e.g., n > scale-free models in our examples) seem to show a small deviation 
from the log-normal prediction at the tails of the distribution, especially for 8 £ — 0.5. In 
particular, the log-normal PDF reproduces very well the skewness and kurtosis measured 
from the simulation data, when the finite size of the simulation volume is properly taken 
into account. 

The degree of agreement of the log-normal models that we have shown is amazing 
considering the fact that the underlying mapping between the initial and the evolved density 
fields differs significantly from the simulation results even in an averaged sense. We have 
explicitly shown the probability distribution of the initial and the evolved density fields from 
simulations, although we were not able to provide a physical explanation for the origin of 
the log-normal PDF. This should be left as our future work and we would like to come back 
later elsewhere. For this purpose, other theoretical approaches based on perturbation theory 
(Bernardeau 1992, 1994) and the spherical collapse model (Fosalba & Gaztahaga 1998) may 
be helpful. 
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Nevertheless our present work provides an empirical justification for the use of the 
log-normal PDF in a variety of theoretical model predictions. For instance, Matsubara & 
Yokoyama (1996) proposed to evaluate the effect of the nonlinear gravitational evolution on 
the genus statistics using the log-normal mapping. Taruya & Suto (2000) constructed an 
analytical model for halo biasing on the basis of the one-point log-normal PDF of underly- 
ing mass density field. Hikage, Taruya, & Suto (2001) applied this biasing model in their 
predictions of the genus for clusters of galaxies. 

Finally, the present results might be useful in considering the prediction of weak lensing 
statistics (Valageas 2000; Munshi & Jain 2000). To construct a model for PDF in redshift 
space is another important topic (e.g., Watts & Taylor 2001; Hui, Kofman & Shandarin 
2000), which is relevant in discussing Ly-a forests (Gaztahaga & Croft 1999). 
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Table 1. Simulation parameters for the CDM models. 



Model 


^0 


Ao 


rt 


0"8 


^boxf/^Mpc] 


realizations 


SCDM 


1.0 


0.0 


0.50 


0.6 


100 


3 


LCDM 


0.3 


0.7 


0.21 


1.0 


100 


3 


LCDM300 


0.3 


0.7 


0.21 


1.0 


300 


3 


OCDM 


0.3 


0.0 


0.25 


1.0 


100 


3 



^Shape parameter of the power spectrum. 



Table 2. Amplitude of a Q \(R) evaluated from the CDM simulations. The values in 
parenthesis denote those estimated from the nonlinear fitting formula of Peacock & Dodds 

(1996). 



smoothing 


R [h^Mpc] 


SCDM 


LCDM 


OCDM 




2 


2.33 (2.24) 


4.17 (4.08) 


4.37 (4.23) 


top-hat 


6 


0.79 (0.77) 


1.37 (1.40) 


1.37 (1.38) 




18 


0.23 (0.24) 


0.44 (0.50) 


0.43 (0.47) 




2 


1.11 (1.08) 


1.95 (1.96) 


1.97 (1.96) 


Gaussian 


6 


0.36 (0.35) 


0.64 (0.69) 


0.63 (0.67) 




18 


0.065 (0.090) 


0.15 (0.22) 


0.13 (0.21) 
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Table 3. Rms a n \(R) in the scale-free simulations. 



smoothing 


R [Lhox\ 


n = 1 


n = 


n = -1 


n = -2 




0.02 


2.48 


3.10 


3.18 


2.79 


top-hat 


0.05 


0.80 


1.10 


1.28 


1.28 




0.15 


0.15 


0.28 


0.43 


0.54 




0.02 


1.00 


1.34 


1.51 


1.47 


Gaussian 


0.05 


0.26 


0.44 


0.62 


0.71 




0.15 


0.03 


0.09 


0.17 


0.25 



Table 4. Amplitude of £, n \(r;R) evaluated from the CDM simulations with Gaussian 
smoothing. The values in parenthesis denote those estimated from the nonlinear fitting 

formula of Peacock & Dodds (1996). 



R [/i^Mpc] 


r [h^Mpc] 


SCDM 


LCDM 


OCDM 


2 


4 


0.68 (0.64) 


2.15 (2.10) 


2.05 (2.06) 


2 


6 


0.36 (0.35) 


1.10 (1.15) 


1.06 (1.10) 


6 


12 


0.058 (0.063) 


0.21 (0.26) 


0.20 (0.24) 


6 


18 


0.021 (0.028) 


0.10 (0.144) 


0.087 (0.127) 



Table 5. Two-point correlation £ n i(r; R) in the scale-free simulations with Gaussian 

smoothing. 



R [Lbax] r [L hox ] n = I n = n = -1 n = -2 



0.02 0.04 0.40 0.82 1.29 1.35 
0.02 0.06 0.13 0.36 0.67 0.83 
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Table 6. Two-point moments, {(SiS 2 ) 2 ) and ((8i5 2 ) s ), from simulations and the log-normal 
PDF predictions. The LCDM models adopt the Gaussian smoothing with R = 2/i~ 1 Mpc 
and the moments are evaluated at the pair-separation of r = 4/i _1 Mpc. For the scale-free 

models, R = 0.02L box and r = 0.04L box . 



Model 


<(W 2 > 


<(W 3 > 


^min 




LCDM 


(0.76 ± 0.4) x 10 3 


(1.5 ± 1.4) x 10 6 






log-normal 


0.64 x 10 3 


1.1 x 10 6 


-0.96 


90 


LCDM300 


(1.4 ±0.2) x 10 3 


(4.9 ±1.5) x 10 6 






log-normal 


1.3 x 10 3 


6.2 x 10 6 


-0.98 


165 


n = 1 


3.1 ±0.1 


73 ± 11 






log-normal 


4.0 


140 


-0.98 


14 


n = 


24 ±3 


(2.6 ±0.6) x 10 3 






log-normal 


30 


5.6 x 10 3 


-0.99 


30 


n — — 1 


88 ± 18 


(2.2 ± 1.0) x 10 4 






log-normal 


102 


3.9 x 10 4 


-0.98 


42 


n = -2 


120 ± 40 


(4.1 ±2.1) x 10 4 






log-normal 


120 


5.3 x 10 4 


-0.96 


46 
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0.1 1 10 0.1 1 10 100 

1+fi 1+5 



Fig. 1. — One-point PDFs in CDM models with Gaussian (left panels) and top-hat (right 
panels) smoothing windows; R = 2h~ 1 Mpc(cyan) , 6/i _1 Mpc(re<i), and 18/?, _1 Mpc(g , reen). 
The top, middle and bottom panels correspond to the PDFs in SCDM, LCDM, and OCDM. 
The solid and long-dashed lines represent the log-normal PDF adopting a n \ calculated directly 
from the simulations and estimated from the nonlinear fitting formula of Peacock Sz Dodds 
(1996), respectively. The values of <r n i in each panel are estimated from the simulations. 
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1 + 5 1 + 5 

Fig. 2. — Same as Fig 1, but in the scale- free models (n = 1 to —2 from top to bottom)] 
R = 0.02L box (cyan), 0.05L box (red), and 0.15L box (green). 
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i i/i i i i i 




0.1 1 0.1 1 10 



(7 , O , 

nl nl 

Fig. 3. — Normalized skewness S and normalized kurtosis K against a n \ from simulations 
and the log-normal PDF predictions. The symbols represent the values estimated from sim- 
ulations (the quoted la error-bars represent the scatter in the realizations). The meaning of 
predictions plotted in different lines is explained in the text. Top-hat smoothing is assumed. 



0.1 1 10 0.1 1 10 100 

l + (5j l + 6 1 




0.1 1 10 0.1 1 10 100 

1 + 5 j 1+<5j 



Fig. 4. — Two-point PDFs in LCDM model with Gaussian smoothing over R = 
2h~ 1 Mpc(upper four panels) and R = 6/i~ 1 Mpc(7ou>er four panels). The results at sepa- 
ration r = 2R and 3R are plotted. The solid and long-dashed lines represent the log-normal 
PDF adopting a n \ and £ n i calculated directly from the simulations and estimated from the 
nonlinear fitting formula of Peacock & Dodds (1996), respectively. 
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Fig. 5. — Two-point PDFs for different models and smoothing window functions. The upper 
four panels plot the results in CDM models, while the lower four panels in scale-free models. 
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Fig. 6. — Two-point PDFs for the scale-free models at negative (left) and positive (right) 
tails of the distribution of 5 2 . 
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-0.5 0.5 -1 1 



(5(z = 9) 6(z=9) 

Fig. 7. — Contour plots of the joint probability P(5(z) \ 5(z = 9)) in the LCDM model 
(I/box = 300/i _1 Mpc) with Gaussian smoothing window; R = 6/i _1 Mpc (right) and R = 
2/i~ 1 Mpc (right). The top, middle and bottom panels correspond to correlations of S(z = 
2.2), S(z = 1), and 5(z = 0) against 5(z = 9), respectively. The white lines represent the 
log-normal transformation (32) and the magenta lines are the conditional mean at a fixed 
5(z = 9). 
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Fig. 8. — Conditional probability P(5(z) \ 5(z = 9)) for a fixed 5(z = 9) corresponding to 
each panel of Fig. 7. In each panel, results for S(z = 9) = —2, —1, 0, 1, 2, and 3 times the 
a(z = 9), the rms of 6(z = 9), are plotted. 



